Scattering of electrons in graphene by clusters of impurities 
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It is shown that formation of clusters of charged impurities on graphene can suppress their con- 
tribution to the resistivity by a factor of the order of the number of impurities per cluster. The 
dependence of conductivity on carrier concentration remains linear. In the regime where the cluster 
size is large in comparison to the Fermi wavelength, the scattering cross section shows sharp res- 
onances as a function of incident angle and electron wavevector. In this regime, due to dominant 
contribution of scattering by small angles, the transport cross section can be much smaller than the 
total one, which may be checked experimentally by comparison of the Dingle temperature to the 
electron mean free path. 

PACS numbers: 72.10.-d; 72.80.Rj; 73.61.Wp 



I. INTRODUCTION. 

Graphene currently attracts intense attention as a 
novel, strictly two-dimensional system with unique elec- 
tronic properties that are interesting with respect to both 
basic physics and potential applications (for review, see 
RefsUj.2.3). It was shown already in the early reports 
on graphene 4 . that charge carriers in this material exhib- 
ited a remarkably high mobility /i such that submicron 
mean free paths were routinely achievable, and an or- 
der of magnitude higher /i were observed for suspended 
graphene samples^). Away from the neutrality point, 
the conductivity of graphene is weakly temperature de- 
pendent and approximately proportional to the carrier 
concentration nZ&. Despite extensive experimental and 
theoretical efforts, there is still no consensus about the 
scattering mechanism limiting fi in graphene on a sub- 
strate. Charged impurities are probably the simplest 
and thus the most natural candidat e) 9 ' 10 ' 11 , and this con- 
jecture is in agreement with the experiments in which 
potassium atoms were deposited on graphene at cryo- 
genic temperatures 1 ^. However, room-temperature ex- 
periments with gaseous adsorbates such as NO2 have 
showed only a weak dependence of ji on charged impurity 
concentration^ 3 -. The latter observation agrees with sev- 
eral reports of only modest changes observed in /i after 
thermal annealing of spuriously doped samples. Further- 
more, recent experiments 1 ^ did not find any significant 
dependence of fj, on immersing graphene devices in high- 
k media such as ethanol and water (dielectric constants 
k ss 25 and 80, respectively) but this also disagrees with 
another report 1 - 5 in which two monolayers of ice increased 
/X in graphene by « 30%. Because of the experimental 
controversy, alternative mechanisms such as scattering 



on frozen ripples 1 ^ and resonant impuritie o 17 ' 18 were dis- 
cussed. 

Regardless of the experimental debate about the dom- 
inant scattering mechanism, the case of graphene cov- 
ered with adsorbates at elevated temperatures^ gener- 
ally requires more careful consideration since there is a 
vast literature which shows the formation of clusters of 
different metals on the surface of graphit e 19 ' 20 ' 21 ' 22 ' 23 ' 24 . 
These atoms easily diffuse on graphite's surface overcom- 
ing only relatively low barriers, and tend to form clusters. 
Potassium atoms on graphite arrange themselves into the 
so called p(2 x 2) structure, with a K-K spacing of 0.492 
nm, that is, roughly, 3.5 nearest-neighbor carbon-carbon 
distances^ 2 -. However, in the case of graphite, this usu- 
ally happens only at low temperatures and high cover- 
age by adsorbates 2 ^. For low doping concentrations such 
as those used in typical experiments on graphene, adsor- 
bates on graphite are randomly dispersed and, at elevated 
temperatures, evaporate from its surface, except for such 
materials as, for example, Au, that forms stable clusters 
on graphite. 

From this surface science perspective, graphene is dif- 
ferent from graphite, and we expect that clusters can be 
more easily formed on graphene and be stable at high 
temperatures. Indeed, it was shown experimentally 13 
that graphene binds such molecules as NO2, NH2, 
H2O, etc. even at room temperature. In the case 
of graphite, they can attach only below liquid-nitrogen 
temperatures 2 ^,. The reason for the stronger attachment 
remains unclear but could be due to the presence of rip- 
ples on graphene 2 ^. According to both experiments and 
theorySS,, ripples can bind even atomic hydrogen that is 
unstable on a flat surface on both graphene and graphite. 

We believe that, once attached to graphene (and this 
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certainly happens for various gases even at room tem- 
perature), adsorbates should tend to cluster, much more 
so than for the case of graphite's surface. First, ripples 
would obviously force them to move from the valleys onto 
the hills which favor the adsorption. Second, there exists 
an additional long range attraction due to Casimir-likc 
interaction mediated by Dirac fermions2£, which is ab- 
sent for graphite. 

On the basis of the above consideration that agrees 
with what is now known about graphene adsorbates, 
both theoretically and experimentally, it is important to 
consider how such clustering of adsorbates can influence 
the electronic properties of graphene. In this report, we 
analyze the scattering of Dirac fermions by clusters of 
charged impurities and show that for the same doping 
level such a disorder results in significantly lower resis- 
tivity. This model reconciles the doping experiments at 
cryogenic*^ and ambient^ conditions, as low tempera- 
tures prevent surface diffusion and, therefore, clustering 
of adsorbates. 

The next section presents the model to be studied. Sec- 
tion III contains the main results. We discuss in section 
IV possible extensions of the model. The main conclu- 
sions are described in Section V. 



II. THE MODEL 

Let us first assume that the charged impurities inside 
the cluster are ordered occupying positions over the cen- 
ters of carbon hexagons, as in the p(2 x 2) structure men- 
tioned above^ 2 .. In such a situation the impurities do not 
break the sublattice symmetry and cannot lead therefore 
to the gap opening. The main effect is therefore merely 
a local doping of graphene, that is, shift of its chemical 
potential, similar to what happens for graphene on the 
top of metals 2 ^. Another effect, that is, the residual un- 
screened Coulomb potential, of the cluster as a whole, 
~ 1/r, far from the cluster, will be discussed further. 

We start with the simplest model, that is, the scat- 
tering of the charge carriers by a closed region where 
the chemical potential has been modified. For simplic- 
ity, we assume a circular cluster. The problem of scat- 
tering of the 2D massless Dirac electrons by the circu- 
larly symmetric potential well has been considered in 
Refs. 17.29.30.31.32.33. The model parameters are the 
Fermi energy and Fermi wavevector outside the cluster, 
ep and kjr, the change in chemical potential inside the 
cluster, V, the Fermi velocity, vp, and the radius of the 
cluster, R. We take ft = 1 in the following. The differen- 
tial cross section can be written in terms of Bessel func- 
tions, whose dimensionless arguments are <f) ou t — kpR 
and fan — {hp + V/vp)R. We assume that the cluster is 
heavily doped, so that 4>i n ^> 4> ou t- The charge induced 
inside the cluster is estimated as n(VR) 2 Vp 2 cx </> 2 „. We 
will neglect intervalley scattering, which is justified if the 
boundaries of the cluster are smooth on the atomic scale, 
and R 3> a, where a is the lattice constant. 




FIG. 1: Angular dependence of the cross section, u(9), in 
nanometers, for a cluster of radius R — 20 nm with a 
chemical potential of V = 500 meV. Red: charge density 
p = 5 x 10 12 cm~ 2 (£; F = 250meV, k F R = 7.9). Green: An- 
gular dependence of the cross section (multiplied by 100) for 
p = 10 10 cm- 2 {E F = llmeV, k F R = 0.35). 



III. RESULTS. 



The scattering cross section read a 17 ' 29 i 30 ' 31 i 32 ' 33 : 
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Note that, since R n = i?_i_ n , the back-scattering am- 
plitude vanishes, f(9 — tt) = which is the consequence 
of the pseudospin conservation at the "crural" scattering 
related with the Klein paradox 34 . 

The cross section shows two regimes, depending on 
whether (j> out — kpR <C 1 or <fi out 3> 1. In the first case, 
the cluster is small compared to the Fermi wavelength. 
The cluster perturbs weakly the electronic wavefunctions, 
and the Born approximation can be used. The differen- 
tial cross section, a (6) has in this case a weak depen- 
dence on the scattering angle 9. The total cross section 
increases as kp is increases, a ~ [V '/ \v pR )] 2 kpR 2 . 

For kpR ^> 1, the cross section as function of the in- 
cident angle 9 shows a narrow maximum at 9 = 0. In 
addition, both the angular resolved and the integrated 
cross sections show resonances, associated to quasi-bound 
states inside the cluster. The integrated cross section de- 
cays slowly as a function of kp. The angular dependence 
of the cross section is shown in Fig. [TJ 

Results for the transport cross section, a t r = 
J 7 ^ 7T a(9)[l — cos(9)]d0, are shown in Fig. [21 We analyze 
in Fig. [5] the total cross section for t V — 0.5 eV, which 
describes the shift in chemical potential due to weakly 
coupled adsorbates, like Al, Ag, or Cr* 2 ^. Similar results, 
although with a smaller periodicity, are found for V = 2 
eV, which describes strongly coupled adsorbates, such as 
K, where the charge transfer can reach 1/8 per carbon 
atom2£. The radius of the cluster was chosen as R = 20 
nm, which is comparable to the size in ripples found in 
graphene 25 . The total number of electrons inside the 
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FIG. 2: Integrated transport cross section cr tr , for a cluster 
of radius R = 20 nm and a shift in the chemical potential of 
V = 0.5 eV. 



high angular momenta in the quantum model. These or- 
bits are probably less efficient in modifying the scattering 
process than the orbits with a lower number of internal 
reflections, leading to the calculated cross section, with a 
sharp maximum as function of the incident angle. Note 
that the resonances under discussion are two-dimensional 
analogs of the "Fabry-Perot" resonances in the Klein tun- 
neling regime!™ 

The elastic electron mean free path, I, is given, approx- 
imately by 



cluster is therefore Ni n = npR 2 « 250, where p is the 
charge density inside the cluster, p = kp /ir,v F kp = V. 

The limit k F R 3> 1 can be analyzed by using the 
asymptotic expressions for the Bessel functions at x — > 
oo: 



J n (x)+iY n (x)*J— e^-^-f) 



(2) 



Then, the expression for the reflection coefficient in radial 
waves, r n (see Ref. [HI) simplifies to 

r w I tan(0;„ - 4> out ) = tan {j^j n < <t>in /g\ 
[ n > <j> in 

and the cross section can be approximated as 
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where n max 
regime is 



fc^-R- The transport cross section in this 
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The scattering process in this limit can be studied by 
the methods of geometrical optic a 32 ' 36 ! 37 . Typical tra- 
jectories, as a function of the shift in potential inside 
the cluster and impact angle are shown in Figs. [3lThc 
scattering will be dominated by periodic orbits inside 
the cluster. These periodic orbits are the semiclassical 
analogues of the resonances of the quantum model. For 
energies such that the internal trajectories are not peri- 
odic, the transmitted waves will interfere destructively. A 
periodic trajectory will lead to transmitted rays at well 
defined angles, as found in the full calculation of a(8). 
Typical trajectories, as function of the shift in potential 
inside the cluster and impact angle are shown in Figs. [3] 
and |U The only periodic orbits for large values of V/E F 
include many internal reflections, which correspond to 



where nc is the cluster concentration. At low carrier 
densities, kpR <C 1 the Born approximation gives: 



<7tr oc kpR 
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and atr is proportional to the density of states and to the 
square of the potential. At high densities, kpR ^> 1, one 
can use Eq.([5]). The conductivity is estimated as 
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(8) 



We expect the oscillations of the cross section shown in 
Fig. [2] to be averaged out in clusters with less symmetric 
shapes. The parameter kpR reaches the value kpR « 
10 — 12 for R = 20 nm and charge density in the clean 
regions p = 2 x 10 13 cm -2 . 

Interestingly, for the regime 4> out 3> 1 the total cross 
section a tot distinguished from <j t r by the absence of the 
factor 1 — cos 9 in Eq. ([5]) is larger than <r tr by a factor 
kpR. The total cross section is related with the single- 
particle decoherence time which determines, e.g., Dingle 
temperature in the Shubnikov - de Haas oscillations 3 ^. 

The elastic mean free path depends on the cluster den- 
sity and carrier concentration, p. For nc — 10 10 cm -2 
and p = 5 x 10 12 cm -2 we obtain I = l/(nco") ~ 200 nm. 

We have neglected so far the long-range part of the 
Coulomb potential induced by the cluster. This poten- 
tial will modify the scattering cross section for electron 
wavelengths kp 1 > R. The cross section will depend on 
carrier concentration as atr oc fc^ 131 ' 39 ' 40 ' 41 . As a result, 
we expect that the conductivity for k F R <C 1 will scale as 
fc 2 ! 11 ' 42 , instead of the dependence given by Eq. ((SJ . How- 
ever, since the scattering cross section is proportional, in 
Born approximation, to the charge square and to the first 
power of the charge concentration, the clusterization will 
lead to suppression of this contribution to the resistivity 
by a factor of order of number of atoms in cluster, in 
comparison with the case of chaotically distributed im- 
purities. 
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FIG. 3: Classical trajectories of an electron scattered by a 
circular cluster. 50 internal reflections are shown. The impact 
angle, 6, of the incoming trajectories is 9 — 7r/4. Left: V = 
E F . Right: V = 10E F . 




FIG. 4: As in Fig. [3] for V = WEf, as function of the impact 
angle. Left: 9 = tt/20. Right: 9 = tt/2 - tt/20. 



IV. 



BEYOND THE SIMPLIFIED MODEL 



Our model of completely ordered impurities inside the 
cluster is oversimplified. However, if disorder inside the 
cluster is relatively weak so that the local mean free 
path I exceeds the electron wavelength inside the cluster 
A w hvp/V one can expect that above consideration is 
correct, at least, qualitatively (the local mean free path 
is defined here as the mean free path of infinite disor- 
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FIG. 5: Angular dependence of the cross section when the 
cluster is determined by a mass term, which breaks the sym- 
metry between the sublattices. The parameters used are 
kR = 10 and AR/v F = 20. 



dered system with the same chemical potential and the 
same distribution of the scattering potential as inside the 
cluster). If the disorder becomes stronger one reaches at 
some moment the Mott limit I « A without further local- 
ization, due to the Klein tunneling^!. In this regime, the 
electron rays inside the cluster are no more straight and 
the "Fabry-Perot" resonances are destroyed. The cluster 
with such strong disorder will behave just as an obstacle 
of size R, with the transport cross section of order of R. 

Another effect which should be considered is a pos- 
sible formation of superstructure inside the cluster. It 
can break the sublatticc equivalence and lead to the lo- 
cal gap opening. To see potential consequences of this 
local reconstruction of the electronic structure one can 
extend the model to the case when the cluster is defined 
by a mass term rather than a shift of chemical potential. 
Similar boundary conditions were discussed in Ref. 0. 
We assume that the mass, A, is only finite inside the 
cluster. We also neglect, for simplicity, by the shift of 
the chemical potential. The cross section in such model 
is expressed in terms of the new reflection amplitudes (cf. 
Eq.©): 



ia-J n (k F R)I n+ i(KR) - a + J n+1 (k F R)I n (KR) 
ia-Y n (k F R)I n+1 (nR) - a + Y n+1 {k F R)I n (nR) 



^A 2 - (E F + Vf 
v F 



2 2E 



(9) 



where I n (x) is a modified Bessel function, which is zero 
at the origin and grows exponentially as x — > oo. 

We have also calculated the cross section including a 
staggered potential, A. The main effect of a mass term 
seems to be to reduce the oscillations of the transport 
cross section as a function of angle. If the mass term is 
large enough, the effect should be qualitatively the same 
as for the strong disorder, that is, the transport cross 
section will be of the order of R, as for a nontransparent 
obstacle in optics. The changes induced by a mass term 
in the differential cross section are shown in Fig. [5] 



V. CONCLUSIONS. 

Let us summarize the main results of our consideration. 

(i) The transport cross section of charge carriers in 
graphene by large neutral clusters due to a shift of the 
chemical potential inside the cluster becomes indepen- 
dent of the cluster size, R, and shift in chemical poten- 
tial, V, for k F R 3> 1, except for an oscillatory func- 
tion. This can be viewed as a consequence of the Klein 
tunneling 3 electrons can always tunnel into the cluster, 
irrespective of the value of V. The oscillatory function 
is, most likely, replaced by its average for clusters with 
irregular shapes, as one can assume by analogy with the 
geometric optics^. 

(ii) The total scattering cross section, obtained by inte- 
grating a (6) over angles, is proportional to R for k F R 3> 
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1, as it should. In this regime a to t/^tr ~ kpR^S> 1 which, 
is principle, can be observed by comparison of the mean 
free path with the Dingle temperature if this scattering 
mechanism is dominant. For all other scattering mech- 
anisms considered before, including charged impurities, 
<Jtot ~ ctr, with a numerical factor of order of one. 

(hi) The transport cross section is proportional to kp , 
Hence, scattering by large clusters leads to a dependence 
on carrier density similar to that for charged impurities 
or resonant scatterers, g oc n. 

(iv) The main difference in the expression for the con- 
ductivity between scattering by neutral clusters and scat- 
tering by charged impurities is that the impurity con- 
centration has to be replaced by the cluster concentra- 
tion which increases the electron mobility, roughly, by 
two orders of magnitude. Thus, possible clusterization 
of charged impurities in graphene can probably explain 
the relatively weak dependence of the mobility on charge 
impurity concentration^ and dielectric constant^. 

v) The formation of clusters is a process favored by 
high atomic diffusion. Hence, we expect that, by anneal- 
ing the samples used ir*i£ above 100K the mobility will 
increase towards the values measured before doping by 
potassium. 

vi) The correlation observed ir>2£ between the shift of 
the Dirac point and the electron mobility for different 
adsorbates as a function of adsorbate concentration is 
consistent with the formation of clusters. The effective 
charge, q*, transferred from the adsorbate atom to the 



graphene layer varies for different adsorbates. For el- 
ements that transfer to graphene an amount of charge 
much less than one electron charge, such as Pt, the scat- 
tering cross section^ goes as q* 2 . The shift of the Dirac 
point should be proportional to oc riiq* , where rii is 
the concentration of the adsorbate. The change in mo- 
bility should scale, on the other hand, as /x" 1 oc riiq* 2 . 
For adsorbates such that q* w le, the mobility scales 
as oc riiq*. Hence, different adsorbates should show 
different ratios [i~ 1 /ED i . A ratio that varies weakly for 
different adsorbates is more consistent with the existence 
of clusters, each of which transfers to graphene a few free 
electron charges, independently of the type of adsorbate 
and size of the cluster. 
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